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ABSTRACT 

If the Pioneer Anomaly (PA) was a genuine dynamical effect of gravitational origin, 
it should also affect the orbital motions of the solar system's bodies moving in the 
space regions in which the PA manifested itself in its presently known form, i.e. as 
a constant and uniform acceleration approximately directed towards the Sun with a 
non-zero magnitude A Pio = (8.74 ± 1.33) x 1(T 10 m s~ 2 after 20 au from the Sun. 
In this paper we preliminarily investigate its effects on the orbital motions of the 
Neptunian satellites Triton, Nereid and Proteus, located at about 30 au from the 
Sun, both analytically and numerically. Extensive observational records covering sev- 
eral orbital revolutions have recently been analyzed for them, notably improving the 
knowledge of their orbits. Both analytical and numerical calculations, limited to the 
direct, Neptune-satellite interaction, show that the peak-to-peak amplitudes of the 
PA-induced radial, transverse and out-of-plane perturbations over one century are 
up to 300 km, 600 km, 8 m for Triton, 17,500 km, 35,000 km, 800 km for Nereid, 
and 60 km, 120 km, 30 m for Proteus. The corresponding orbital uncertainties ob- 
tained from a recent analysis of all the data available for the satellites considered 
are, in general, smaller by one-two orders of magnitude, although obtained without 
modeling a Pioneer-like extra-force. Further investigations based on a re-processing 
of the satellites' real or simulated data with modified equations of motions including 
an additional Pioneer-type force as well are worth being implemented and may shed 
further light on this important issue. 
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and satellites: individual (Neptune, Nereid, Proteus, Triton) 



INTRODUCTION 



The Pioneer Anomaly (PA) jNietol [2006) consists of an 
unmodeled, almost constant and uniform acceleration ap- 
proximately directed toward s the Sun and of magnitude 
jAnderson et al.lll99Sl , l2002aT ) 



A 



.74 ± 1.33) x 10~ 10 m s" 2 . 



(1) 



It was detected in the radiometric data from the Pioneer 
10/11 spacecraft after they passed the 20 au threshold mov- 
ing along roughly antiparallel escape hyperbolic paths taken 
after their previous encounters with Jupiter (~ 5 au) and 
Saturn (~ 10 au), respectively. The PA's existence has been 
subsequently confi r med by inde pend ent investigations by 
iMarkwardtl (|2002l ), lOlser] (|2007h and iLevv et"aTT (|2009T ) as 
well. Interestingly, latest data analyses are focussing on pe- 
riodic variations of the anomaly, characterized as functions 
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of the azimuthal angle tp defined by the direction s Sun- 
Earth antenna and Sun-Pioneer l|Levv et al.l I2009D . Con- 
cerning the possibility that it s tarted to manifest itself 
at shorter heliocentric distances (.Ni eto fc Anderson 120051 ; 
lNietdl2008l '). efforts to retrieve and analy ze early data from 
Pioneer 10/11 are current ly being made (|Toth fc Turvshevl 
|200S| ; iList fc Mullinl lioosh . The PA is one of some astro- 
metric anomalies in the solar s ystem reported in recent 
years (Lammerzahl et all 120081 ; lAnderson fc Nietol l20ld : 
Iloridl2009al ~ 

Attempts to explain some features of the PA in terms 
of mundane, non-gravitational effects, pertaining the Pio- 
neer probes themselves lik e thermal forces among differ - 
ent p arts of the s pacecraft jAnderson et aL 2002a; iMurph 



1999: Katz 1999 



iScheffer] [^OoiTTMbelek 



Michalskil 1200: 

Bertolami et all 120081 : iToth fc Turvshevl I2009T) or externa l 
influences like anisotropic solar emission (jBini et all 120041 ), 
have been undertaken, but some of them have n ot obtained 
full consensus so far l|Anderson et all Il999allbl: I Biui et al 
l2004h . On the other hand, latest work by iBertolami et al 
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(2008) strongly points out that PA is a thermal effect due 
to the energy sources in the spacecraft; further studies on 
possible thermal effects like a potential asym metric heat dis- 
sipat ion of the spacecraft surface are ongoing (|Rievers at al.l 
2009). Conventional explanations of gravitational origin in 
terms of drag due to interplanetary dust, dark matter , 



Kuiper Belt Objects ([Anderson et al 



Nieto et alj|2005t Ide Diego et~ 



200 



2002al; iNietol 120051; IRathke fe Izzol (|200< 



By fitting the TNOs' observatio ns with modi fi ed eq uations 
of motion according to eq. |T]), Wallin et all (|2007f ) found 
(0.87 ± 1.6) x 10~ 10 m s~ 2 , which is consistent with zero 
and whose upper bo u nd is i nconsistent with e q. (QT ) at 4cr 
level. IRathke fc Izzol d200d) , Ilorio fc Giudicel l|20oi ). Ilorid 

looked at the outer planets. 



2007al) 



and 



Tangcn 



Bertolami fc Vieiral 



200d ) have been found not satisfactorily as well. As a con- 



sequence, many suggestions invoking non-standard gravita- 
tional and non-gravitational p hysics like non-linear electro- 
dynamics (Mb elek et al.ll2007l) have been proposed. For a 
review see, e.g., (lAnderson et all |2002al; IDittus et all 120051: 
Bertolami fc Paramos! 120061 : IRathke fc Izzol 1200(1 " de Diegj 

ous pro- 



20081 ) and references therein. Among the 



pos ed exotic gravitation a l me chanisms we recall those 
by iBrownstein fc MoffatJ (|2006l ). based on a long-range 
Yukawa-like extra-force, and by iJaekel fc Revnaudl (|2008l ) 
who proposed metric extensions of the Einstein's Gen- 
eral Theory of Relativity (GTR). Attempts to find ex- 
otic gravitational explanations for PA did not even cease 
after the publication of the latest works on the non- 
gravitational effects like lBertolami et al.l (|2008 | ): just to lim - 



iting to published works , see, e.g 



IWilson fc Blomel ||2009|) :Ig reaves 



Avramidi fc Fuccil 
(2009) 



Exirifard 



2009) 



2009) 



A dedicated spaceraft-based mission to test the PA in the 
outer regions of t he solar system has also been proposed 



outer regions oi th e solar s yst em nas also been proposed 
and investigated (IDittus et all 120051 : IRathke fc Izicl I2OO6I ; 
iBertolami fc Pa ramos 2007). 



The hypothesis that non-standard forces of gravita- 
tional origin are able to explain the anomalous behavior of 
the Pioneer spacecraft must cope with the following cru- 
cial remark. If the PA was due to some modifications of the 
known laws of gravity, this should be due to a radial extra- 
force affecting the orbits of the astronomical bodies (planets 
and their satellites, comets, Trans-Neptunian Objects, etc.) 
as well, especially those moving in the space regions in which 
the PA manifested itself in its presently known form. Oth- 
erwise, a violation of the equivalence principle largely in- 
compatible with the present-day boun ds of ~ 10~ 13 from 
Earth -based laboratory experiments ( Sch lamminger et all 
2008) would occur. The impact of a Pioneer-like additional 
acceleration on the motion of planets and minor bodies 
in the outer regions of the solar system interested by the 
PA was recentl y studied by numerous authors with differ- 
ent approaches (lAnderson et al.ll2002bl; I Rathke fc IzzollgOOfj; 



Iorio fc Giudicel l2006l; I Page et al.ll2006l: iPitieval l2006j ; Ilorid 
2007allbl; iTangenl 120071; IWallin et all 120071; Standisbl l200gf 
Ioridl2009d; IPage et al.l 120091: [Fienga et al.ll2009l ). In par- 



ticular, lAnderson et all (2002b) discussed the impact of a 



Pioneer-like acceleration on the long-period comets and the 
form of the Oort cloud; however, such bodies are not par- 
ticularly well suited to perform accurate tests of gravita- 
tional theories because of the impact of several aliasing 
non-gravitation al perturbations l ike out-gassing as they ap- 
proach the Sun. lPage et all |2006l ) investigated the potential 
offered by an analysis of the minor planets in the outer solar 
system to con firm or refute the e xistence of a gravitational 
Pioneer effect. IWallin et all (|2007i ) used a well-observed sam- 
ple of Trans-Neptunian Objects (TNOs) between 20 and 
100 au from the Sun to constrain Pioneer-like deviations 
from Newtonian gravity in that region of the solar system. 



parameterized the PA in terms of a 
change of the effective reduced solar mass felt by Neptune 
finding it nearly two orders of magnitude beyond the current 
observational constraint. Moreover, they noted that the Pi- 
oneer 11 data contradict the Uranus ephemerides-obtained 
without explicitly modeling the PA-by more than on e or- 
der of magnitu de. Ilorio fc Giudicel l|2006l ). Ilorid l|2007al ) and 
ITanged (2007) computed the secular effects induced by an 
unform and radial extra-acceleration like that of eq. {TJ on 
the orbits of Uranus, Neptune and Pluto, located at 20-40 
au from the Sun, and co mpared them to the p resen t-day, 
unmod ified ephemerides. Ilorio fc Giudicd (|2006D and Ilorid 
(2007a|) concluded that the resulting anomalous effects on 
all of them would be too large to have escaped from de- 
tection sojai. Doubts concerning Neptune were raised by 
ITanged (|2007l ) in the sense that the accuracy of the cur- 
rently available observations for it would not, in fact, ex- 
clude the possibility that Neptune is acted upon by Api a . 
Other authors made a step further by including a Pioneer- 
like extra-acceleration in the force models and fitting again 
the planetary observations with such modifie d equations of 
motion. More specifically, IPage et all (|2009l ) fitted modi- 
fied dynamical models including eq. {TJ to observational 
records for Uranus, Neptune and Pluto showing that the 
current ephemeris of Pluto does not preclude the existence 
of the Pioneer effect because its orbit would not be well 
enoug h charac t erized at present to make such an asser- 
tion. [Standis^ (|2008l ) fitted planetary data records with a 
modified version of the JPL DE ephemerides with a uni- 
form extra-acceleration directed towards the Sun acting on 
Uranus, Neptune and Pluto; a magnitude as small as just 
10% of eq. {TJ yielded compl etely unaccep t able r esiduals for 
all the three outer planets. iFienga et all (|2009l ) added an 
extra-acceleration like that of eq. {1} to the equations of mo- 
tion of the outer planets and fitted the resulting modified 
ephemerides to their observations by finding that Uranus 
excludes the existence of Pioneer-like acceleration as large 
as eq. JTJ) at a 4<r level. On the contrary, for Neptune and 
Pluto the effect of eq. ([TJ is absorbed by the fit, so that 
the resulting residuals do not allow to exclude the existence 
of a Pioneer-like anomalous acceleration affecting such bod- 
ies. The existence of a standard PA in the r egions crossed 
by J upiter an d Satu rn has been ruled out by Ilorid |2007d ) 
and IStandishl (|2008h with different approaches. For non- 
standard, velocity-dependent forms of the PA and their com- 
patibility with differe nt ep hemer i des of t he outer planets, see 
IStandishl <| 20081 , 1 20 id ) and Ilorid (|2009bl ); such different ap- 
proaches show that almost all of them are not compatible 
with the planetary observations. 

In this paper we investigate a different astronomical lab- 
oratory with respect to those examined so far to put on 
the test the hypothesis of the gravitational origin of the 
PA independently of what detected in the Pioneer 10/11 
telemetry. Indeed, we will look at the orbital effects in- 
duced by a Pioneer-like acceleration directed towards the 
Sun on the Neptunian satellites Triton, Nereid and Proteus 
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Table 1. Semimajor axis a, eccentricity e, inclination /, longi- 
tude of the ascending node f!, argument of pericenter lo and or- 
bital period f\, of Triton, Nereid and Proteus. Reference frame: 
ICRF/J2000.0 with Neptune as center body. Reference system: 
Ecliptic and Mean Equinox of Reference Epoch. Epoch: 31 Octo- 
ber 1989. The HORIZON interface by NASA has been used. The 
orbital period of Neptune is 164.9 yr. 



Keplerian orbital element 


Triton 


Nereid 


Proteus 


a (km) 


354, 767 


5,517,147 


117,714 


e 


0.00003 


0.75428 


0.00090 


/ (deg) 


130.9 


5.0 


28.9 


a (deg) 


213.2 


320.3 


48.1 


w (deg) 


60.2 


296.1 


54.4 


Pt (d) 


5.8 


360.4 


1.1 



Table 2. Cartesian coordinates and distance of the Sun, in au, at 
the epoch. Reference frame: ICRF/J2000.0 with Neptune as cen- 
ter body. Reference system: Ecliptic and Mean Equinox of Ref- 
erence Epoch. Epoch: 31 October 1989. The HORIZON interface 
by NASA has been used. 



20 (au) 


1/0 (au) 


z e (au) 


r Q (au) 


-6.2 


29.5 


-0.4 


30.1 



in vi ew of the recent improvements in their orbit determina- 
tion (|Jacobson| :2009) based on the analysis of extensive data 
records covering several orbital revolutions. Their Keplerian 
orbital elements are listed in Table Q] In this paper we will 
perform a preliminary sensitivity analysis by means of ana- 
lytical and numerical calculations. Their goal is to check if 
the scenario considered is worth further, more detailed in- 
vestigations. They could involve, e.g., a re-processing of the 
Neptunian satellites' real or simulated data sets with mod- 
ified equations of motion including a standard gravitational 
Pioneer-like extra-acceleration radially directed towards the 
Sun as well. 

In Section [2] we first analytically work out the anoma- 
lous PA-type orbital effects on Triton, Nereid and Proteus 
(Section 12. Then, we perform numerical integrations of 
their equations of motion with and without the PA. Finally, 
we compare our results to the latest determinations of the 
orbital accuracies for such satellites (Section 12. 2) . Section [3] 
is devoted to the conclusions. 



2 EFFECTS OF A STANDARD PIONEER 
ANOMALY ON THE NEPTUNE'S 
SATELLITES 

We use ICRF/J2000.0 with Nept une as cen t er bod y as refer- 
ence frame. To be consistent with ljacobsonl (|2009) , we adopt 
31 October 1989 as reference epoch. The reference system 
used has the ecliptic and mean equinox of reference epoch. 
In such a frame a standard Pioneer-like acceleration has the 
form 



A p, ° = A 1 



n , 



(2) 



where uq is the unit vector pointing towards the Sun dis- 
played in Table [2] As a consequence, A Pl ° has the compo- 



Table 3. Components of A pio = A pio n = 
A Plo /ro{:E0, j/0, 2q}, in m s~ 2 , at the epoch. Reference 
frame: ICRF/J2000.0 with Neptune as center body. Reference 
system: Ecliptic and Mean Equinox of Reference Epoch. Epoch: 
31 October 1989. 

A Pi ° (m s" 2 ) A Pio (m s" 2 ) A Pi ° (m s~ 2 ) 



-1.79 x 10- 



8.55 x 10" 10 



-0.13 x 10- 



nents shown in Table [3] at the reference epoch. It can be 
noted that it is mainly directed along the y axis of the cho- 
sen frame. Since we are interested in its secular, i.e. aver- 
aged over one orbital period, effects on the motion of the 
Neptunian satellites, we can safely consider A Pl ° as con- 
stant because of the short satellites' periods (see Table [1} 
with respect to the Neptunian one amounting to 164.9 yr. 
In other words, each satellite faces the action of a constant 
and uniform disturbing acceleration directed along a generic 
direction in space which, in general, does not coincide with 
the Neptune-satellite radial one. 



2.1 Analytical and numerical calculation 

The orbital effects of such an anomalous acceleration can be 
worked out with standard perturbative techniques by using, 
e.g., the Gauss equa tions for the variations of the elements 
|Bertotti et alj|2003h 



dt 

de 
dt 

dl 
~dt 

dn 

dt 
If 



— [eAflSin f + A T (- 
nn L V r 



V_ 

na 



A R sin / + A T 



cos / + 



— — An (—') cos it, 
nan V a / 

An (J- J sin it, 



na sin In 



V 
nae 



-Ar cos / + At 1 + - sin / 



dM 2 . ( r \ ( dui T dQ 

— : — = n Ar I — — r) — — hcosi— — 

dt na \a / \ dt dt 



(3) 
(4) 
(5) 
(6) 

coslH^) 
(8) 



where M is the mean anomaly of the orbit of the test 
particle, / is its true anomaly reckoned from the peri- 
centre position, u = w + / is the argument of latitude, 
n = sj GM/a :i — 2ir/Pb is the unperturbed Keplerian 
mean motion (G is the Newtonian constant of gravitation 
and M is the mass of the central body), r\ = \/l — e 2 and 
p = a(l — e 2 ) is the semi-latus rectum. Ar, At, An are the 
projections of the perturbing acceleration A onto the radial 
R, transverse T and out-of-plane directions of the parti- 
cle's co-moving frame whose time-varying unit vectors are 
l|Montenbruck and Gillll2000l ) 



(cos Q cos u — cos / sin Q, sin u 
sin Q cos u + cos / cos Q, sin u 
sin / sin u 



(9) 
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(— sin it cos 0. — cos I sin Q cos u \ 
— sin sin u + cos 7 cos Q cos it (10) 
sin I cos u y 

(sin I sin f2 \ 
— sin I cos O (11) 
cos I J 

A straightforward calculation shows that the R — T — N 
components of a constant and uniform perturbing acceler- 
ation, like our A Pl ° over the timescales involved here, are 
linear combinations of A x , A y , A z with coefficients propor- 
tional to harmonic functions whose arguments are, in turn, 
linear combinations of u, Q and /.Thus, in the satellite's 
co-moving frame A Pl ° is time-dependent through / in u. 
In order to have the secular perturbations of the Keplerian 
orbital elements, Ah, At, An have to be inserted into the 
right-hand-sides of the Gauss equations which must be eval- 
uated onto the unperturbed Keplerian ellipse 



r — a(l — e cos E), 



(12) 



where E is the eccentric anomaly, and integrated over a full 
orbital revolution by means of 



~K 

Other useful relations are 

cos E — e 

cos f — — , 

J 1 - e cos E ' 



(1 — ecos E) 



2tt 



sin / = 



dE. 



Vl-e 2 sin£ 
1 — e cos E 



(13) 



(14) 



After cumbersome calculations it turns out that, apart from 
the semimajor axis a whose secular rate vanishes, all the 
other Keplerian orbital elements ip experience non-vanishing 
secular precessions of the fornfl 



fW a 4. r yw> A 4- C KW> A 



In it. 



ip = I,tt,u,M. (15) 



7 MO 



x,y,z, 



(16) 



where F^\e) are complicated functions of the eccentric- 
ity and are linear combinations of the longitude of the 
pericenter w = oj + f2, f2 and I. In principle, they are time- 
varying according to 



w = vjq + zut, 
Q = Q + fit, 
/ = Jo + It; 



(17) 
(18) 
(19) 



from a practical point of view, since their secular rates are 
quite smaller, especially for Triton and Nereid, we can as- 
sume tjj m zuq, Q. w Oo, / ~ Io in computing cos £jk , where 
vj, Qo, Io are their values at epoch (see Table [IJ. 

The R — T — N shifts over a gen eric time interv al At 
can be exactly computed according to ICasottol |l993). The 
radial perturbation is 



AR = K a Aa + K e Ae + Km AM, 



(20) 



where 
K a = 

Ke = 

K M = 



-a cos /, 
ae 



sin/. 



(21) 
(22) 
(23) 



eq. (|20[) shows that it would be incorrect to identify the shift 
in the radial component of the orbit with the perturbation 
of the semimajor axis only. Otherwise, misleading conclu- 
sions concerning the mean motion n and, thus, the trans- 
verse component as well could be traced. Indeed, if, say, a 
secular signature in AR was found, from the identification 
AR = Aa it could be argued that an analogous perturbation 
in the mean motion 

^=-5^£ (24) 
n 2 a y J 

would occur as well. As a consequence, a quadratic effect 
in the transverse component should occur through the per- 
turbed mean longitude. Actually, this does not happen in 
our case: indeed, we will show that, although, as already 
noted, no secular effects on a are present, both the radial 
and the transverse components exhibit cumulative pertur- 
bations with secular trends, without any quadratic signature 
in the transverse one. The transverse perturbation is 



AT = H e Ae + HmAM + r(Aw + cos/Afi), 



with 



1 + 



sin /, 



H 



M 



a 2 VT 



(25) 

(26) 
(27) 



The out-of-plane perturbation is 

AN = r (AI sin u — Af2 sin I cos u) . (28) 

In eq. (|20[l . eq. (|25jl and eq. (|28p the perturbations of the 
Keplerian orbital elements have to be intended as 

Ai>= f dip, ip = a,e,I,Q.,uj,M, (29) 
Jo 

where dip is taken from eq. (|3}-eq. @. The explicit ex- 
pressions of Ai? pio , AT pi °, AiV pi ° are rather cumbersome, 
so that we will not explicitly show them. It turns out that 
linearly growing signatures are present in all of them along 
with sinusoidal terms; quadratic terms are, instead, absent. 

In Figure [l]-Figure [3] we plot the Pioneer- induced R — 
T — N perturbations for Triton, Nereid and Proteus over 
a century. In order to express the eccentric anomaly as a 
function of time we used a partial sum of the seriefl 



E = M + 2J2 



Js(se) 



sin(sM), M = n(t — t ), (30) 



where J 3 (se) are the Bessel functions of the first kind. 

The peak-to-peak amplitudes for the Pioneer-type R — 
T — N perturbations are 300 km, 600 km, 8 m for Triton, 
17, 500 km, 35, 000 km, 800 km for Nereid, and 60 km, 120 
km, 30 m for Proteus. 



1 The semimajor axis a does not appear in the denominator of 
the equation for the eccentricity rate. 



2 It converges for all e < 1. See on the 
http://mathworld.wolfram.com/KcplersEquation.html 
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Triton 



Nereid 




Figure 1. Analytically computed transverse, radial and out-of- 
plane Pioneer-induced shifts for Triton over 100 yr according to 
eq. (|206 . eq. l|25|l and eq. {28) , and Table [T] In the expansion of 
eq. (1306 we retained just the first term because, in view of the 
extremely small eccentricity of the orbit of Triton, the second 
term is of the order of 10 — 10 . 
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Figure 2. Analytically computed transverse, radial and out-of- 
plane Pioneer-induced shifts for Nereid over 100 yr according to 
eq. \20\ . eq. j25j> and eq. {28) , and Table [T] In the expansion of 
eq. J 30 1) we retained the first 40 terms because, in view of the 
large eccentricity of the orbit of Triton, the following ones are of 
the order of, or smaller than 10 -6 . 



It should, now, be pointed out that we have only con- 
sidered the direct perturbations induced by A Pl ° on each 
satellite considered separately; in fact, the total effect may 
be even larger because of the mutual gravitational interac- 
tions among the satellites themselves and Uranus which are 
all allegedly influenced by the PA as well. Moreover, it can be 
argued that, over time intervals larger than one orbital pe- 
riod as those used here, the Pioneer-induced signatures are 
modulated by the slowly changing Q,u,I because of Nep- 
tune's oblateness and N— body interactions with the other 
giant planets and satellites themselves, and by the variation 



of the Sun's position which reflects into slow changes in the 
components of A Pl °. 

To support our analytical calculation and to further 
clarify the issue of the semimajor axis, we also performed nu- 
merical integrations with MATHEMATICA of the equations 
of motion of Triton, Nereid and Proteus with and without 
the PA. Concerning the classical forces common to all the 
three satellites, we included the first two even zonal harmon- 
ics J2, J4 of Neptune and the attraction of Uranus, Saturn 
and Jupiter. The intersatellite interactions have been taken 
into account as well by considering Nereid and Proteus as 
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Proteus 




-40 -20 20 40 60 

t(y) 
Proteus 




-40 -20 20 40 60 

t(y) 

Proteus 
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Figure 3. Analytically computed transverse, radial and out-of- 
plane Pioneer-induced shifts for Proteus over 100 yr according to 
eq. J20J, cq. l|25|l and eq. {28) , and Table [T] In the expansion of 
eq. (1306 we retained just the first term because, in view of the 
extremely small eccentricity of the orbit of Triton, the second 
term is of the order of 10 — 7 . 



massless point par ti cles a cted upon by a massive Triton, as 
done by iJacobsonl (2009). First, in Figure [4] we plot the 
Pioneer-type perturbations on the semimajor axes a of Tri- 
ton, Nereid and Proteus over one Keplerian orbital period. 
As expected from our analytical calculation, no cumulative, 
net effects occur; this is neither in contradiction with the 
absence of a quadratic signature in AT nor with the pres- 
ence of a linear signature in AR. Then, for each satellite we 
computed Ar(t) = r pi °(t)-r Ncvton (t) and projected it onto 
the R — T — N co-moving frame of the unperturbed orbit 
to have AR pio (t), AT pio (i), AN pio (t). It turns out that our 



numerical integrations confirm the analytical calculations: 
for saving space we do not show here the pictures of the 
numerical integrations. 

2.2 Confrontation with the accuracy of the orbits 

Although obtained differently, Fig ure [T}Figure E] h ave can 
be compared with Figures 4-6 by IJacobsonl (|2009r ). which 
yield a measure of the accuracy of the orbits of the Neptu- 
nian satellites considered. They have been obtained by fit- 
ting the observational data sets without including the PA in 
the dynamical force models of Triton, Nereid and Proteus, 
but are useful to give us an idea about a possible detec- 
tion of Pioneer-type effect s in their orbital dynamics. More 
specifically, Figures 4-6 bv lJacobsonl (|2009h display the orbit 
uncertainties from the consider covariance mapped into the 
R — T — N directions. In general, the solution covariance 
yields an optimistic measure of the orbit uncertainties be- 
cause it does not account for possible systematic or unmod- 
eled errors. They can occur in the dynamical force models, in 
the obser v ation modeling or in the observation themselves; 
IJacobsonl l|2009h believes that the dominant systematic er- 
rors mainly reside in the observations because the models 
adopted fit to them at their pr esumed ac c uraci es. To include 
the effect of neglected errors, IJacobsonl (|2009T l added some 
"consider" parameters to the estimation process. They are 
quantities which are not estimated, but whose uncertainty 
contributes to the uncertainty in the estimated parameters. 
Th e reliability of su ch a procedure was subsequently tested 
by IJacobsonl £2009) refitting the orbits with different data 
sets and comparing the consequent changes in the orbits to 
the uncertainties derived from the consider covariance for 
the modified data sets. It turned out that the changes in 
the orbits were at or below the level of uncertainties. Thus, 
IJacobsonl l|2009r ) concluded that all important errors were 
properly accounted for, and that the consider covariance can 
reliably be adopted as a realistic measure of the orbit accu- 
racy. 

Figure 4 by IJacobsonl l|2009h deals with Triton. The 
middle panel shows the radial distance uncertainty which 
oscillates between 0.8 km and 1.2 km over a time span of 
one century. The middle panel of Figure Q] tells us that 
the Pioneer-induced radial perturbation would be as large 
as 300 km for Triton, so that it seems reasonable to argue 
that the resulting overall anomalous shift should not have 
escaped from detection over a time span of more than one 
century; note that the astrometry of Triton covers 161 yr 
from 1847, i.e. one year after its discovery, through 2008. 
Even by re-scaling the radial uncertainty by a factor 10, 
the situ ation wou l d not change. The upper panel of Fig- 
ure 4 bv lJacobso n (2009) depicts the transverse uncertainty. 
It linearly grows reaching a level of about 150 km after 60 
yr; the upper panel of Figure [1] shows that AT Pl ° is larger, 
although not by two orders of magnitude as in the radial 
case. Thus, also a Pioneer-type transverse effect may have 
remained u ndetectable with difficulty. The lower panel of 
Figure 4 bv lJacobsonl (120091 )) displays the out-of-plane un- 
certainty which amounts to about 50 km after 60 yr, while 
the anomalous PA perturbation would be orders of magni- 
tude smaller, as shown by the lower panel of Figure [T] 

The o rbit a ccuracy of Nereid is shown in Figure 5 by 
IJacobsonl (|2009l ). The radial distance uncertainty linearly 
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grows up t o about 1 , 600 k m after 60 yr (middle panel of 
Figure 5 bv lJacobsonl ( |2009h ). while the peak-to-peak radial 
PA effect is as large as Ai? Pl ° = 17, 500 km for Nereid whose 
data set covers 59 yr from its discovery in 1949 through 2008, 
i.e. one o rder of m a gnitud e larger. The upper panel of Fig- 
ure 5 by I Jacobsonl l|2009l ) displays the transverse accuracy 
which linearly grows up to about 3, 000 km. According to 
Figure [2l AT Pl ° = 35, 000 km, i.e. about one order of mag- 
nitude larger. The accuracy in the out-of- plane direction i s 
displayed in the lower panel of Figure 5 bv lJacobsonl (2009); 
it oscillates between a few km to 60 km. According to the 
lower panel of Figure [2] , the corresponding Pioneer-type 
out-of-plane peak-to-peak amplitude is as large as 800 km, 
more than ten times larger. Thus, in the case of Nereid the 
PA perturbations over 59 yr would be more than one order 
of magnitude larger than the corresponding orbit accuracy 
in all the three directions. 

Figure 6 bv lJacobsonl (|2009T) depicts the T - R - N ac- 
curacy for Proteus whose observational record is about 20 
yr long. The radial dis tance uncertainty , shown in the mid- 
dle panel Figure 6 by Ijacobsonl ([2009) , oscillates between 
2.5 km and 4.5 km, while the peak-to-peak amplitude of 
the anomalous PA radial signal over a similar time span is 
about 20 — 30 km. The tran sverse accuracy (upper panel of 
Figure 6 by Ijacobsonl (|200gh ) linearly grows up to 200 km, 
making, thus, problematic a detection of the corresponding 
PA transverse shift (AT Pl ° = 120 km). The anomalous out- 
of-plane effect is several orders of magnitude smaller than 
the correspon ding orbit accura cy shown in the lower panel 
of Figure 6 bv lJacobsonl (|2009l ') which is as large as 50 km. 
Thus, in the case of Proteus a PA-type radial perturbation 
would be about one order of magnitude larger than the cor- 
responding orbit uncertainty, while the transverse and out- 
of-plane PA signatures would have been overwhelmed by the 
corresponding orbit uncertainties. 

Finally, let us conclude by noting that, concerning the 
direct e ffect of the PA on the orbital motions of the outer 
planets. Ilorio &: Giudicd l|2006l ) showed in their Table 1 that 
the induced anomalous perihelion precessions xbpio are in the 
range 83.5 — 116.2 arcsec cty -1 for Uranus-Pluto. Actually, 
latest determinations of t he corrections Aw to the standard 
perihleion precessions bv lPitieval l|2010l) with the EPM2008 
ephemerides summarized in her Table 8 are, instead, — 3.89± 
3.90 arcsec cty -1 , -4.44 ± 5.40 arcsec cty -1 , 2.84 ± 4.51 
arcsec cty -1 , respectively. 



3 DISCUSSION AND CONCLUSIONS 

We have investigated the impact that an anomalous, con- 
stant and uniform acceleration directed towards the Sun 
having the same magnitude of the PA would have on the 
orbital dynamics of the Neptunian satellites Triton, Nereid 
and Proteus which move in the deep PA region of the solar 
system. Long data sets covering a large number of orbital 
revolutions are currently available for them. 

We, first, used an analytical approach which only con- 
sidered the direct PA-type perturbations on the three satel- 
lites taken separately to work out the corresponding shifts 
in the radial, transverse and out-of-plane orbit components. 
In fact, also the indirect effects caused by the PA-affected 
mutual gravitational interactions among them should be, in 



principle, considered. Then, we numerically integrated the 
equations of motion with and without an extra-PA accel- 
eration confirming the analytical findings. It turned out 
that only secular and sinusoidal signatures are present in 
the three orbit components; we showed that this is not in 
contrast with the absence of a secular effect on the semi- 
major axis. No quadratic terms appear in the transverse 
component, as, instead, it would happen if a was affected 
by secular signatures. 

Our analysis showed that the resulting anomalous or- 
bital effects are much larger than the realistic orbit accu- 
racies evaluated from a recent analysis of all the available 
astrometric observations by one-two orders of magnitude. 
However, it must be stressed that our investigation should 
be considered preliminary. Indeed, it would be necessary 
to refit the entire set of observations to the corresponding 
predictions computed by taking the anomalous PA effect 
into account. As an alternative approach, it would also be 
possible to fit the predicted observations without the PA to 
a set of simulated observations produced by including the 
PA. Our study demonstrates that such further investiga- 
tions, which are beyond the scopes of this paper, should be 
considered worth the needed effort. 
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Figure 4. Pioneer-induced shifts on the semimajor axes for Tri- 
ton, Nereid and Proteus. They have been obtained by numer- 
ically integrating with MATHEMATICA the equations of mo- 
tions of the three satellites with and without the PA in the 
form j4 Plo Ti0 starting from the same set of initial conditions at 
the epoch 31 October 1989. Then, we have computed Aa(t) = 
a pi °(t) - a Ncwton (t). The gravitational attractions by Uranus, 
Saturn and Jupiter and the effect of the first two even zonals 
of the non-spherical gravitational field of Neptune have been in- 
cluded. Given the relatively short time interval of the integration, 
the positions of Uranus, Saturn, Jupiter and the Sun have been 
kept fixed to their values at the epoch computed with the HORI- 
ZONS software. 



